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Abstract 

We consider a class of one-dimensional Heisenberg spin models (plaquette 
chains) related to the recently found 1/5-depleted square-lattice Heisenberg 
system CaV^Og. A number of exact and exact-numerical results concerning 
the properties of the competing dimer and resonating plaquette phases are 
presented. 
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Recent discovery of a quantum disordered phase and a spin gap in the quasi-two- 
dimensional magnet CaV40Q has attracted considerable interest The magnetic 
system can be described by the antiferromagnetic Heisenberg model on a 1/5-depleted square 
lattice containing interacting plaquette spins and additional frustrating interactions ||^. In 
the limit of weak interplaquette and frustrating couplings, the ground state of the model is a 
product of resonating four-spin singlet states formed on each of the elementary plaquettes of 
the lattice. In the other hmit of strong interplaquette couplings, the ground state approaches 
the spin-Peierls phase constructed of dimers connecting neighboring plaquettes. A number 
of interesting questions has arised concerning the nature of the transition between these two 
phases, the possible existence of an intermediate Neel phase, and the role of geometrical 
frustration in the formation of a spin gap. In this paper we address similar problems in 
plaquette Heisenberg chains. Fig. 1, which are closely related to the parent two-dimensional 
depleted lattice described above. Apart from the pure theoretical interest, we believe that 
the plaquette spin chains may appear as real objects in some magnetic systems. To the best 
of our knowledge, plaquette Heisenberg chains have not been studied except for the recent 
work of Katoh and Imada who presented a second-order perturbation estimate for the gap 
of the unfrustrated model in the limit of weak interplaquette interaction We consider 
the simplest Heisenberg plaquette chain described by the Hamiltonian 

H = H, + H^ + H, = J,Y. SiSj + J2 E SiSj + J3 E SiSj, (1) 

where the sums run over the dimer Ji, plaquette J2, and diagonal J3 bonds, respectively. 
Fig. la. For the sake of clarity, we neglect the second diagonal coupling in the four-spin 
plaquettes. The exchange constants Ji, J25 and J3 are supposed to be non-negative. The 
so-defined model is perhaps among the simplest ones displaying the tendency to four-spin 
plaquette and spin-Peierls dimer ground states in the limits of weak and strong interplaquette 
(dimer) interaction Ji, respectively. The inclusion of H-^ in Eq. (1) opens a chance to study 
the role of frustration in the formation of of a spin gap structure and the relative stability 
of these two phases. 
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Note also that the square of the the local spin operator Soi = Sai + S;,i ( i = 1, 2, Np, 
Np = N/ 4) is a good quantum number, so that the " diagonal" plaquette spin Soi defined by 

= Soi{soi + 1) can be used to classify the states of the chain by the set [soi, So2, SoNp] 
defining a sector of states for every given Soi '■ \ Sai — Sbi |< Soi < Sai + Sbi- Here Sai and Sbi 
are the spins of the diagonal sites a and b in the elementary plaquette i. Due to the relation 

+ 1) = 

+ 1) + Sbi{Sbi + 1) + 2SaiSbi, (2) 

the local operator SaiSbi (and the global one H^), introducing a geometrical frustration in the 
plaquette i, is diagonal and explicitly known in every subspace connected to [soi, So2, SoNp]- 
In the particular case Sai = Sbi = s = 1/2, the number of sectors corresponds to the 
number of all possible sets [11 10 10.. 01], where 1(0) stands for Soi = l(soi = 0). 
(i) Unfrustrated plaquette chain, J3 = 0: 

An exact statement concerning the structure of the ground state of the unfrustrated 
J3 = plaquette chain can be deduced from the Lieb-Mattis theorem for Heisenberg 
antiferromagnets on bipartite lattices. A corollary of the latter theorem [0, valid for the 
ground state from the subspace with a given total spin, reads (SjSj) > provided that the 
pair of sites (i,j) belongs to one and the same sublattice, and (SjSj) < in the opposite 
case. Applied to the diagonal spins in a given plaquette, the sign rule reads {SaiSbi) > 0. 
Thus, for a s = 1/2 plaquette chain one obtains the following statement. The ground state 
of a s = 1/2 plaquette chain without frustration is a singlet state and belongs to the sector 
[11. ..1]. In the general case of arbitrary site spins, the above sign rule, combined with Eq. 
(2), implies 

Soi{Soi + l) > Sai{Sai + l) + Sbi{Sbi + l), i = l,...,Np. (3) 

One can check that for Sai = Sbi = s and s > 2, the diagonal plaquette spin Soi has more 
than one possible value for every plaquette. 

Now, let us consider unfrustrated s = 1/2 open chains containing Np plaquettes and two 
end spins of the form shown in Fig. la. We already know that the ground state belongs 
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to the sector [11...1]. There is a simple way to build classes of singlet excited states using 
the following construction. Consider the lowest excited state of the chain in the sector with 
exactly one Soi = plaquette, say, [110111] for Np = 6. Due to the presence of a singlet 
dimer bond in one of the plaquettes, the state is a product of the ground states of two 
independent clusters of the same form ( containing Np = 2 and Np = 3 plaquettes in our 
example). If the Soi = plaquette is situated in the end of the chain, then the product will 
contain the state of a A'p — 1 plaquette cluster and two singlet dimer states formed on the 
(ab) and (cd) end bonds. This factorization feature reflects the fact that the third spin in 
a triangle is completely free provided that the other two spins are in a singlet dimer state. 
We can continue the dividing procedure by adding new Soi = plaquettes. The final state, 
belonging to the sector [000000], is a spin-Peierls dimer state. Note that the above classes 
of states are directly connected to the states of the periodic chain, namely, the ground state 
of a periodic chain with Np plaquettes in the subspace with one soi = plaquette is just the 
ground state of an open cluster with A'p — 1 plaquettes times the singlet dimer state on the 
i-th (ab) bond. The above construction is applicable without changes also to the frustrated 
system. 

In Fig. 2 we present exact diagonalization (ED) results (for periodic chains) for the 
correlations of the nearest-neighbor spins on the plaquette (ac) and dimer (cd) bonds. It is 
seen that already at J2 = Ji the resonating plaquette character of the short-range order is 
practically completely established. In the region near J2 ~ 0.5Ji, the competition between 
the dimer and the plaquette type short-range correlations is strongly pronounced. Notice, 
however, that the symmetry of the ground state is not changed for arbitrary finite J2. 

In Fig. 4 are in particular shown ED data for the spin gap function A(J2) of the 
unfrustrated model (open circles). For small enough J21 the numerical data predict A ~ J|, 
whereas for J2 > 1 one has A ~ J2- Here Ji = 1. In the small J 2 limit, the above behavior 
can also be deduced from the perturbation theory. Alternatevely, for small enough J2 one 
can also apply the well-known Dasgupta-Ma eliminating procedure [Q in order to remove 
the dimer spins from the Hamiltonian. The resulting effective model, on the other hand, was 
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shown by Schulz to be equivalent to the 5" = 1 Heisenberg chain. Thus, the finite spin- 
gap structure in a s = 1/2 plaquette chain for small J2 corresponds to a picture where the 
low-energy physics of the model is dominated by the dynamics of the "solid" diagonal spins 
Soi coupled by an effective antiferromagnetic exchange interaction Jg// ~ Jl/J^ produced 
by the spins on the dimer bonds, i.e., an effective Haldane-type model with a spin gap 
A ~ Je//exp(— vrso). So = 1 [0]. Close to the dimer limit J2/J1 ^ 1, the ED data show 
spin-spin correlations characteristic of a spin-one Heisenberg antiferromagnetic chain. 
{ii) Frustrated plaquette chain, J3 > 0: 

Below we proceed to frustrated periodic chains containing Np elementary plaquettes. It 
is easy to show that for sufficiently large J3 the dimer eigenstate 

1*0 >= 11x1? (4) 

i=l 

is an exact ground state of the model. Here xi?(0 = Wa{i)l3h{i) — l3a{i)oib{i)]l and 
Xcdi^) = [o^cii — ~ l3c{i — l)ttd(0]/v^5 <^ (3 being the usual = ±1/2 states for 

s = 1/2. This can be achieved most readily if one prove that the energy of the above state 

E,iNp) = ~J^Np-p3Np (5) 

saturates a lower bound to the ground state energy for suitable values of the parameters Ji, 
J2, and J3. To this end, one can use the simplest cluster decomposition of the Hamiltonian 
H = Yl'i=i hi, where hi = (Ji/2)SiS2 + ^282(83 + S4) + (^3/2)8384 is the Hamiltonian of 
the cluster shown in Fig. lb. The variational principle implies the following relations for 
the ground state energy Eo{Np) 

EdiNp) > EoiNp) > 2Npe, (6) 

where e = -(Ji/8) - (J2/4) + (J3/8) - y^Jf - 2J1J2 -|- 971/4 is the ground state energy of 
the cluster. The latter inequality implies the following bound above which |\E^o > is an exact 
ground state 
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J3c = -2 + f + 2Vl-2'/2 + 9J|, Ji^l. (7) 

In Fig. 3 we compare this bound with ED data for a cluster with N = 24 sites (the data 
for = 16 are very close to the presented numerical data). It is seen a good theoretical 
estimate in view of the fact that we have used the simplest available cluster. 

For the model under consideration one can also use the following expression for the 
ground state energy in the subspace with n diagonal spins in a soi = 1 state and Np — n 
diagonal spins in a soj = state 

^^(iVp) = + Jgn - h^Np, (8) 

where Eq^\Np) is the ground state energy of the unfrustrated Hamiltonian Hq = Hi + H2 
in the same subspace. From the condition E^'^\Np) = Ed{Np) one can find the coupling J3 
at which the energy levels ^(")(A/p) cross Ed{Np). In the case n = Np one gets the following 
bound 

MN,) = 'I ~ J. - 1. (9) 

where Eq{Np) is the ground state energy of the J3 = periodic chain with Np plaquettes. 
In Fig. 3 is shown that the latter expression for a periodic chain with Np = 6 plaquettes 
exactly represents the numerical ED data for the exact bound above which the state |\E'o > 
is a ground state of the model. However, it is not clear if the above expression also describes 
the exact absolute bound in the thermodynamic limit 

In Fig. 4 we show ED data for the spin gap function A(J2) for J3 = 0.5. One can clearly 
see different behavior of A(J2) before and after the point ~ 0.53. The latter point is 
close to the estimate from Fig. 3 at J^^ = 0.5. Notice also that the data in Fig. 4 in the 
region J2 > J2C coincide with the results in the unfrustrated system. The latter observation 
is a result of Eq. (8) provided that the lowest excited state belongs to the sector of the 
ground state. Finally, the ED data show that in the region J2 < J2C the J3 = 0.5 numerical 
data are in essence independent of N for the studied clusters, N = 8, 16, and 24. 
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(Hi) Conclusions: 

The above results demonstrate the role of frustration in the plaquette chain model. In the 
frustrated plaquette chain one finds two different (by symmetry) phases. For large enough 
J2 and fixed J3, the ground state belongs to the sector [11... 1] and it has features which are 
typical of the unfrustrated system. For small J2 plaquette coupling, the ground state belongs 
to the sector [00... 0] which is characterized by a spin-Peierls type dimer state. In principle, 
one should expect a quantum first-order phase transition between the above phases deserving 
a special consideration. In the two-dimensional analog of the plaquette chain, there is a new 
role of the frustration. The latter may also destroy the possible intermediate Neel phase 
which does not present in the one-dimensional system. 
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FIGURES 

FIG. 1. a) The plaquette chain model studied in the paper. Ji, J2, and J3 are respectively 
the dimer, plaquette, and diagonal antiferromagnctic couplings. The spins leaving on ab and cd 
sites in every four-spin plaquette are called diagonal and dimer spins, respectively, b) The cluster 
used to find a lower bound above which the dimer state is an exact ground state of the plaquette 
chain. 

FIG. 2. Exact diagonalization data for the short-range plaquette, < Sa(i)Sc(i) >, and dimer, 
< Sc(i)Sd(i + 1) >, spin-spin correlators respectively normalized by the factors —1/2 and —3/4, 
vs J2. Ji = 1. 

FIG. 3. The bounds Jsc vs J2 above which the state 1^*0 >, Eq. (4), is an exact ground state 
of the frustrated plaquette chain. The solid line represents the analytic result, Eq. (9), for Np = 6. 
The dashed line is the lower bound, Eq. (7), obtained with the cluster shown in Fig. lb. The open 
circles are numerical ED results for a periodic chain with Np = 6 plaquettes, i.e., AT = 24 sites. 
Ji = l. 

FIG. 4. Exact diagonalization results for the spin gap function A(J2) of the J3 = and J3 = 0.5 
plaquette chain. The extrapolation is deduced from the N = 8, 16, and 24 site clusters. Ji = 1. 
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Fig. 1 
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